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Geometric (Clifford) algebra provides an efficient mathematical language for describing physical 
problems. We formulate general relativity in this language. The resulting formalism combines the 
efficiency of differential forms with the straightforwardness of coordinate methods. We focus our 
attention on orthonormal frames and the associated connection bivector, using them to find the 
Schwarzschild and Kerr solutions, along with a detailed exposition of the Petrov types for the Weyl 
tensor. 
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| /" ( ' Geometric (or Clifford) algebra provides a simple and natural language for describing geometric concepts, a point 
which has been argued persuasively by Hestenes Q and Lounesto Q among many others. Geometric algebra (GA) 
unifies many other mathematical formalisms describing specific aspects of geometry, including complex variables, 
, matrix algebra, projective geometry, and differential geometry. Gravitation, which is usually viewed as a geometric 
■ theory, is a natural candidate for translation into the language of geometric algebra. This has been done for some 
aspects of gravitational theory; notably, Hestenes and Sobczyk have shown how geometric algebra greatly simplifies 
certain calculations involving the curvature tensor and provides techniques for classifying the Weyl tensor 0, Q. 
Lasenby, Doran, and Gull [jj have also discussed gravitation using geometric algebra via a reformulation in terms of 
a gauge principle. 

• In this paper, we formulate standard general relativity in terms of geometric algebra. A comprehensive overview 
like the one presented here has not previously appeared in the literature, although unpublished works of Hestenes and 
of Doran take significant steps in this direction. We obtain several useful results not available in standard treatments, 
including a general solution for the connection in terms of coordinate frame vectors (Section [HJ and simplified forms 



I . for the Petrov classification of the Weyl tensor (Section IVIj l. Various mathematical objects, such as the connection 
^jjy and the curvature tensor, gain explicit geometrical interpretations through their mathematical forms. Also, geometric 
algebra provides notable compactness and clarity for many derivations, such as the generalized Schwarzschild solution 
in Section lHI Bl and the Kerr solution in Section lYlII An explicit translation between geometric algebra and differential 
forms is provided in Section ITVl which provides orientation to those unfamiliar with geometric algebra but conversant 
5_j ■ in conventional differential geometry. 

We assume knowledge of Clifford algebra techniques and notation, and their geometric interpretation. Appendix lAl 
reviews basic elements and lists some pedagogical references, while Appendix [B] defines basic derivative operators and 
tensors. Definitions and results from these Appendices will be used without comment. 



II. COVARIANT DERIVATIVES AND CURVATURE 



A. Orthonormal Frames and The Bivector Connection 



Since we are working within the framework of standard general relativity, we will assume the existence of the 
spacetime manifold, and the presence of a metric on it. This means we can define the covariant derivative V in 
the standard way (see Section 3.1 of 0), so that it (1) is linear, (2) obeys the Leibniz rule, (3) commutes with 
contraction, and (4) reduces to d in its action on scalar functions (see Appendix iBl for details); we also impose the 
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torsion-free condition, which we will discuss in more detail later. We specify the operator uniquely by letting it be 
metric-compatible. 

The scalar- valued derivative operator t ■ V can be decomposed into the form 

a- VA = a- VA + r(a,A), (2.1) 

where V is called a connection, which has the same grade (or grades) as A. The symbol V is another derivative 
operator that we can choose: for example, if we want T to be the Christoffel connection, we let V = d. We will 
instead adapt the covariant derivative to a particular frame of vectors, whose "metric" is fixed from point to point in 
the manifold. 

A frame of vectors {e^} is one for which the pseudoscalar (volume element) does not vanish, 

e = e A ei A e 2 A e 3 ^ 0, (2.2) 

for a region of spacetime; the existence of a pseudoscalar in some region is equivalent to the existence of a metric on 
the same region of a manifold Q. (See also Eq. (|A25|> in Appendix El) A given vector in the reciprocal frame is 
found via 

e , = [\ e v e~\ (2.3) 



which satisfies the orthonormality condition 



K- (2-4) 



Our starting point in developing a geometric algebra description of spacetime is a particular class of frame vectors, 
the orthonormal frame {e M } satisfying the condition 

^.e u = r,i" v , (2.5) 

where [77] = diag(l, — 1, — 1, — 1) is the metric for Minkowski spacetime. (In fact, we do not need the frame to be 
orthonormal — all we need is for the metric in this frame to be constant.) We attach this orthonormal frame to each 
point of the spacetime manifold, but we will leave the discussion of metrics and coordinate frames to Section III Dl 
If we apply the covariant derivative to Eq. 12.5fl we find the special property 

a ■ V(e M • e„) = a ■ c?(e M • e v ) = 0; (2.6) 

applying the Leibnitz rule over the inner product yields 

a ■ V(e M • e v ) = (a ■ Ve M ) • e„ + e M ■ (a ■ Ve„) = T a (e M , e v ) + T a (e v , e M ) = (2.7) 

where T is a representation of the connection components with respect to the basis {e^}. Since the function 
r a (e A ,,e ;y ) = — r a (e t/ ,e M ) is antisymmetric and linear in its arguments, we can write it in terms of a bivector as 
follows (as shown in Chapter 3 of Q): 

r a (e M , e„) = n(o) • (e M A e v ) = (11(a) • e M ) • e„, (2.8) 

where fi(a) is the bivector connection. This means the covariant derivative of a frame vector can be written as 

o-Ve^ = n(o).e^ (2.9) 

while for an arbitrary multivector we have the decomposition 

a ■ VA = a ■ VA + fi(a) x A, (2.10) 

where we have imposed the condition a ■ Ve M = for frame vectors and a ■ V</> = a ■ d(f> for scalars; when the 
frame vectors e M are constant, V = d. (Note that x is the commutator product, not the vector cross-product: 
A x B = (AB — BA)/2.) Due to the "fixed" nature of the orthonormal frame, the connection bivector acts as a 
"rotation" of objects with respect to a "background" defined by theie^,}; this is the foundation of the flat-spacetime 
gauge theory of gravity expounded by Lasenby and his colleagues 0, |7j ■ The relation to gauge theory is given in a 
note following this section. 
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With the covariant derivative in hand, we define parallel transport in the usual way; the magnitude of the vector 
being transported is constant, so 

t • V6 2 = = 2b (t • V6) -► t • V6 = 0. (2.11) 

A geodesic is defined when the tangent vector is parallel transported along itself: 

t-Vt = 0. (2.12) 

We can write three derivative operators based on Eq. I|2.10|l : 

V • A = d a ■ (a ■ VA) = V • A + d a ■ (ft (a) x A) covariant divergence (2.13a) 

V A A = d a A (a • VA) = V A A + d a A (O(a) x A) covariant curl (2.13b) 

VA = V ■ A + V A A = d a (a ■ VA) covariant gradient, (2.13c) 

where d a is defined in Appendix [BJ As we will see in Section II VI the covariant curl is equivalent to the exterior 
derivative of differential geometry (when torsion is zero). 



B. Geometric Interpretation of the Connection 



Comparison of two infinitesimally close points in curved spacetime is not as simple as its analogue in a flat space, 
but if we assume that the orthonormal frame may be carried smoothly from one point to another in the manifold, 
many of the results of flat spacetime carry over. The same multivector field at different points A(r) and A(r') along a 
curve f parameterized by some parameter t will not in general have the same magnitude or direction. The covariant 
derivative is a measure of this difference, which can be viewed geometrically as the difference between the multivector 
at a given point on a curve and its value at an infinitesimally displaced point along the curve which has been parallel 
transported (that is, transported without changing its magnitude or direction) back to the original point. Without 
loss of generality we evaluate the derivative at t = 0: 



t ■ VA = lim ~ (/^ r) (A T ) - A T=0 ) = A^ T) (A r ) 



t^O T 



d 

d7-(-^ 



(2.14) 



where 



t[f] 



df(r) 



is a vector tangent to the curve and A,_ T JA T ) is the parallel-transported multivector evaluated at r = 0. We are 

only able to define A in either a flat space, or if we are able to carry the orthonormal frame smoothly along the 
curve f; see 0, Chapter Vbis for more information. 

We can simplify Eq. 12.1411 by letting A7=o = L (the identity operator) and A(r = 0) = A: 



t- VA 



dA d/T 



dr 



dr 



"(A) 



(2.15) 



The type of transformation that does not affect the "absolute value" of a multivector is a (proper, orthochronous) 
Lorentz transformation, which we write as 



A' = A{A) = UAU where UU = 1 
so that A' 2 = A 2 for any multivector A. U is an even multivector called a rotor, which we write as 



U 



3 B/2 



B 



B. 



so that B is a bivector. If we let B(r = 0) = 0, we find that 



(2.16) 
(2.17) 



d/T 



dr 



"(A) 



dU~ TTA dU 
— AU + UA— 

dr dr 



dB 

cb 7 



-A A- 



dB 

d7 



dB 

d7 



x A, 



(2.18) 
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where "x" represents the commutator product (AB — BA)/2. Thus, if we let 

x A ee £l(t) x A, 



dA 

VA= — 

dr 



T = 



dB 

and r(t, A) = — — 
dr 



(2.19) 



r=0 



we have the covariant derivative in the form of Eq. (|2.1|) , with fl being the same bivector connection as in Eq. I|2.10|l . 
The connection bivector maps a vector, which is associated with a path in spacetime, onto a bivector, which is related 
to the (four-dimensional) rotation of the orthonormal frame as it is carried from one point in the manifold to another. 



C. Curvature Tensors and Einstein's Equations 



In general, covariant derivatives will not commute: 



where 



a ■ V, b ■ V A = (a • Vb • V - b ■ Va ■ V) A - [a, b] ■ VA 



[a, b] ee a ■ V6 - b ■ Va 



(2.20) 



(2.21) 



denotes the Lie bracket (we are assuming throughout that torsion vanishes). The accent marks in the first expression 
show that the derivatives are to act only on the multivector, not on the vectors a and b. Using Eq. I|2.10|l for the 
covariant derivatives, we define the Riemann curvature tensor R as 



a ■ V,6- V 



R(a A b) x A, 



R(aAb) ee a ■ VJ7(6) - b ■ Vil(o) + Sl(a) x $7(6) - ft ([a, b]) . 



(2.22) 



This object is a bivector-valued function of bivectors which is linear in both of its arguments, due to the linearity 
of the derivative operators. Also, the outer product guarantees that both sides vanish if a and b are not linearly 
independent. Geometrically, the grade of the curvature tensor and its argument are sensible: the argument specifics a 
particular spacetime plane, and if a multivector is parallel-transported around an infinitesimal closed path contained 
in this plane, it will undergo a rotation described by the bivector returned by the Riemann tensor. 
The Riemann tensor obeys the Ricci identity (see 0, section 11 for derivations) 



the symmetry condition 



and the Bianchi identity 



d a A R(a A b) = 

a ■ R(b A c) + b ■ R(c A a) + c ■ R(a A b) = 0, 



(a A b) ■ R(c Ad) = R{a A b) ■ (c A d), 



V A R(a A b) = 

a ■ V/?(6 A c) + b ■ V/?(c A a) + c • V/?(a A b) = 0. 



(2.23a) 
(2.23b) 



(2.24) 



(2.25a) 
(2.25b) 



In both the Ricci and Bianchi identities, the "a" and "b" expressions are equivalent, although Eq. (|2.23bll and 
Eq. i|2.25b(l are the simpler ones to use in calculations. 

Contraction of the Riemann tensor yields the Ricci tensor: 



R{b) = d a -R(aAb), 
which is vector-valued and has the symmetry properties 

d a A R{a) = 

and 



(2.26) 



(2.27) 



a ■ R(b) = R(a) ■ b. 



(2.28) 
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One last contraction yields the Ricci scalar: 

R=d a R{a). 

The Einstein tensor is defined in terms of the Ricci tensor and scalar as 

G(a) = R{a) - ^aR. 



(2.29) 



(2.30) 



The Bianchi identity, Eq. I|2.25a|) . can then be rewritten in an illustrative form by contracting out the vectors a and 
b: 



d a - 



O b ■ [V ' R(a , b)) = R(V) - ~VR = 



or, more compactly, 



G(V) = 0. 



Finally, the Weyl tensor C(fl A b) is defined by 



C{a /\b) = R(aAb)- - (R(a) A b + a A /?(&)) l^aA bR. 

C is antisymmetric in its arguments, and tractionless: 

C(a A 6) = -C(6Ao) d a C(af\b) = 0. 



Section rVIl is devoted to properties of this tensor. 
The Einstein field equations are 



G(a) = 8irT(a) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



in units with G = 1, where vector- valued T(a) is the stress-energy tensor. For a perfect fluid with density p and 
pressure p, the stress-energy tensor is 



J(EM) ( ) = ( p + p^ a . uu _ p a 



(2.36) 



where p is the density, P is the pressure, and u 2 = 1 is the 4- velocity of an observer comoving with the fluid. Another 
useful form is for electromagnetic fields with field strength bivector F; the stress-energy tensor is then 0, 



V ; 2 



(2.37) 



D. The Coordinate Frame 



After some work, we find the curl of a reciprocal frame vector is 

VAe p = i(e« A e^) ([e a , ep] ■ e») , (2.38) 

which is the Maurer-Cartan equation in curved spacetime 40]. A particularly useful type of frame is the coordinate 
frame {g^} defined (in the absence of torsion) by the holonomic condition 

[fl M , 9 J =0 or <7 M • Vg, = g v ■ Vg^, (2.39) 

which leads immediately to 

VA/E9„A(a. VgV) =VAg" + 9 a A (12(a) • g^ 1 ) - 0. (2.40) 

This Cartan's first structural equation in a disguised form; see Section llVI for an explicit translation of the differential 
form language. Unless otherwise noted, vector and multivector components will be written in terms of the coordinate 
frame; for example, 

b»=b-g^ , V M = g M ■ V, etc. 
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As shown in Appendix [U] we can invert Eq. I|2.40[l to obtain the connection: of a coordinate 

M(b) = ^3"AvV-9 m AW + s"a(6.vs m )) 
n a = n(g a ) = ^{g u AV(g a - g v ) +g» AV a9fi y, 

where V M = g„ • V; the second expression is for the specific case of a frame vector, which is often the more useful 
form. We may write the Riemann tensor over frame vectors as 

R^ = R(g^ A g v ) = g^ ■ V«„ - g v ■ Vft M + ft M x (2.42) 

while the Ricci tensor is 

% = = 9 a ■ R(a A 9/ J = <f • R V(l . (2.43) 
The usual (scalar) components of these tensors are, respectively, 

R\nv = {a K a g x ) ■ Rn» R^v = sv ■ R v ; 

the same process can be used to obtain components of any tensor object, such as the stress-energy tensor of a perfect 
fluid: 

T^ F) = g, • T( PF \g u ) = {p + P)u,u„ - Pg M ■ g v , 

where it„ = u ■ g v . 

Finally, we can assign a linear operator to map the orthonormal frame into the coordinate frame: 

g ll = h- 1 (e tl ) = h;e a . (2.44) 

This operator is a representation of the vierbein field; it occupies a prominent position in the gauge theory formulation 
of gravity given by Lasenby, Doran, and Gull 0, where it acts as a frame map, lifting a frame of vectors in flat spacetime 
toagenerally-covariantspacetimei. The coordinate frame also defines the metric components via 

g^-9u = 9^- (2-45) 

We emphasize that in the formulation outlined here, the metric can be seen as a by-product of the choice of some 
particular coordinate frame. All of the geometrical content of gravitation is defined without reference to the spacetime 
metric g^. In order to define the notion of parallel transport which preserves lengths of vectors, we have alterna- 
tively assumed the existence of an orthonormal frame of vectors at each point in the manifold, in contrast to most 
textbook treatments of general relativity, in which the metric is the central mathematical object (despite its lack of 
direct physical content). However, it should be pointed out that both starting points give mathematically equivalent 
descriptions. 

III. A METHOD FOR DIRECT CALCULATION AND THE SCHWARZSCHILD SOLUTION 

A. General Solution Technique 

Given the constructions of the previous section, we have an explicit expression for the connection bivectors Jl„ in 
terms of the coordinate basis vectors and the curvature bivectors R^ in terms of the connection; thus the Einstein 
equations can be viewed as differential equations in the coordinate frame vectors. Symmetries of the system can be 
incorporated directly into the definition of the coordinate vector. 

For a given physical system which has well-defined symmetries, 

1. Define an orthonormal frame {e^} and its reciprocal {e^} that corresponds to the symmetry at hand (e.g., 
spherical coordinates for spherical symmetry, axisymmetric coordinates if a preferred axis exists, and so on) 

2. Define a frame g^ = /j -1 (e /i ) that preserves the symmetry of item^ 

3. Define the derivative to preserve the orthonormal frame: a ■ Ve M = 0, and decompose the covariant derivative 
as a ■ VA = a ■ V + tt(a) x A; 



frame vector pl|: 

(2.41a) 
(2.41b) 
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4. Calculate the connection bivectors, using Eq. Q2.41b| 

n a = n{g a ) = — (jj v A Vg av + g» A V Q g p 

5. Calculate the Riemann curvature bivectors: 

R»u = ■ - g v ■ vo M + fif, x n u ] 

6. Find the Ricci vectors: 

Rn = g ■ Rvn', 

7. Solve the vacuum equations for the coordinate frame vectors {g v }: 

R„ = 0. 

or the field equations in the presence of some source of energy-momentum: 

1 



where T — ■ T(g ). 



R„ = 8tt I T(g M ) - -g^T 



We propose this technique as a replacement for the usual "brute force" component method, as well as providing an 
alternative to the method of differential forms in cases of special symmetries. As we will see in our discussion of the 
Kerr solution in Section IVlII we often do not even need to specify a particular coordinate system to extract a great 
deal of information from the connection and curvature tensors. This method combines the efficiency of forms with 
the straightforwardness of coordinate methods, and as such we believe can fill the role of either. 

B. Example: Schwarzschild Geometry 

We now use the recipe laid out in the previous section to obtain the exact solution to Einstein's vacuum equations 
for the case of stationary, spherically-symmetric spacetimes. The resulting black hole solution is of course equivalent 
to Schwarzshild's solution, but we provide a more general solution, which in a single expression contains a large class 
of equivalent metric forms, of which the Schwarzshild metric is a special choice. In Section IVIII we will examine the 
more complicated case of the rotating black hole, after we examine techniques to classify and study spacetimes. 

1. Conditions for Spherical Symmetry 

Flat spacetime is trivially invariant under spatial rotations; in fact, every point may be considered the center of a 
sphere for the purposes of rotations. For the non-trivial vacuum case in curved spacetime, we will pick a single point 
as an origin in the orthonormal frame and write a set of basis vectors in spherical coordinates: 

e. il -e u =f\ llv (3.1) 

where = e M (0, <fi) in general, for the angular variables 8 and <fi. We define a rotation as 

R(A) = c^ B Ae^ B (3.2) 

where B = B(x) is a spacelike bivector. The timclike direction is invariant under rotations: B • eo = 0. Our "spherical 
coordinates" reflect this condition in the orthonormal frame. 

The vierbein fields should preserve this spherical symmetry, while allowing for distortion in the radial and temporal 
directions. In other words, we want the rotation operator to commute with the vierbein operation: 

IT 1 (RA) = RhT 1 (A). (3.3) 
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Therefore, define a set of orthogonal basis vectors in spherical coordinates = dx^ for the orthonormal frame: 

e o _ ^,0 e 1 = sin^cos^ 1 + sinflsin07 2 + cos^7 3 ,„ 

e 2 = (cos^cos^ 1 + cos^sin07 2 — sin^7 3 ) e 3 = (r sin#) _1 (— sin 07 1 + cos</>7 2 ) ' 

where the {7^} are the ordinary Cartesian basis vectors for four dimensions, with (7 ) 2 = 1 and (7*) 2 = —1 for 
i = 1,2,3. We normalize these to make their inner product constant: 

e 2 = re 2 and e 3 = rsin#e 3 , (3.5) 

while e° = e° and e 1 = e 1 . Now e M • e v = rf" 1 " ', where 77 = diag(l, —1, —1, —1) is the metric. 

The e° and e 1 directions are invariant under any rotation, while e 2 and e 3 are rotated into each other; any set of 
basis vectors of this type are spherically-symmetric. In particular, a set of basis vectors in the gauge-covariant space 
are automatically spherically-symmetric if they take the following form [{|: 



g* = h(e ) = fre + fa 1 g r = ^(e 1 ) = 9l e x + g 2 e° 

g e = h(e 2 ) = r- 1 (ae 2 + fie 3 ) g* = h(e 3 ) = (r sin^-^ae 3 - f3e 2 " 



9' 


= h(e°) 


= fie + he 1 


9t 


= if 




= .9i e o - 32ei 


g' 


= hie 1 ) 


= qie 1 + g 2 e° 


9r 


= h' 




= fi&i - f 2 e 


9° 


= h(e 2 ) 


-1»2 

= r e 


9o 


= h~ 


-\e 2 ) 
-\e 3 ) 


= re 2 


9 


= h(e 3 ) 


= (r sin^) _1 e 3 


9, t , 


= h~ 


= r sin 8e 3 



where the gi, a, and (3 are functions of r and t. Note, however, that there is no loss in generality in letting f3 = 
and a = 1, since we can effectively perform this operation via rotations and judicious choice of the radial variable r. 

The stationary condition amounts ultimately to the independence of all quantities from the time variable t |(| . Note 
that this is not the same as the static condition, which is a statement that g l ■ g r = 0. The static condition is highly 
restrictive, and leads automatically to the Schwarzschild metric. The stationary condition is sufficient to derive a 
general expression for spherically-symmetric spacetimes (though we need not even assume that — see e.g. 0]), as we 
now show. 

2. Solution of the Field Equations 
Applying the conditions a = 1 and /3 = to Eq. Ij3.6|l yields the following co- and contravariant basis vectors: 



(3.7) 



where the functions fi and gi are assumed to depend on r only, and we have imposed a simple gauge condition 

fx9x - h92 = 1. (3.8) 

This gauge condition is actually the only one we need to solve the vacuum equations; it guarantees that the frames 
in Eq. (|3.7() satisfy the torsion- free condition in Eq. (|2.40|l . 
The directional derivatives of the basis vectors are 

9t ' Y 9 " = g<t> - Vg v = 

9 r ■ Y^t = Sl e ~ 52 e l 9 r ■ = e 2 /g gs 

9 r ■ Vfl r = /i e i - /2 e o flV • ^94, = sin 6e 3 
9e ■ Vg = r cos 8e 3 g g ■ Vg^ = 

where the primes indicate differentiation with respect to r. Plugging these into Eq. (|2.41b|) . we obtain the connection 
bivectors (after some work): 

°t = ~ {gi9i -52ff 2 ) e oei = ~^d r {gl - gl)e e x 

°r = (5l/2 - 92/0 e e l 

tie = gieie 2 - g 2 e e 2 

$7^ = sin 6* (sieie 3 - g 2 e e 3 ) + cos0e 2 e 3 (3.10) 
where we have used the fact that e M A e„ — e M e„ when /i ^ v. Note that no derivatives of the /, functions appear. 
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The six Riemann curvature bivectors follow from Eq. i|2.42[) and take the compact forms 

-^(gf- gi)g t Ag r 
1 a i 2 2^9e^9t 



Rrt 


= g r ■ Vth = 


Ret 


— fig X fl t = 


D 

K<t>t 


= Si0 X fi t -- 


Rer 


= -flr-vn, 


Rfir 


= -g r 'Vn, 


R<pe 


= -g«-vn. 



2 ' Wi r 

1 fl /■ 2 2xflVVSt 

+ fi 9 xfi r = id r ( 5 ?- 52 2 )^^ 



1 2 r 
= -g fl • Vn + 0^xO fl = (( 5 2 _ ^) _ i) (3.11) 

Note that all of these equations contain the quantity g\ — g\ (which is the goo term of the metric). The Ricci curvature 
tensor (via Eq. I|2.43(l ) is particularly simple: 

R t = <f>(r)( gi e - g 2 e 1 ) = <P(r)g t 
R r = $(r)(/iei - / 2 e ) = $(r)g r 
R e = *(r)e 2 - (*(r)/r) flfl 

R* = ^(r)sm0e 3 = (*(r)/r)g (3.12) 

where 

*(') - \^M-9l) + \^9l-9t) = ^Wi-9t)) (3.13a) 
*(r) - - 9l) + \{9l ~ 9l 1) = i {I (r( 5l 2 - <? 2 2 )) - l} • (3.13b) 

In the vacuum case, — identically reducing in our case to $(r) = and \&(r) = 0, which are differential equations 
in the quantity g\ — g\. Letting x( r ) = 9x ~ g\i the differential equations become 

_(r X (r))=0 and — (r X (r)) - 1 = (3.14) 

with the solution 

X(r) = 5? - 5 2 2 = 1 - - (3-15) 
r 

where if is a constant of integration. From past knowledge of the Schwarzschild metric, we know that K = 2M. 

Note that the solution in Eq. 13.15fl and the gauge condition from Eq. I|3.8[) are the only restrictions the vacuum 
equations can give us. We cannot specify from the preceding analysis alone what any of the quantities fi and g{ are 
by themselves. Thus, the most general line element satisfying the field equations is 

dr 2 - 9fl ■ g v dx»dx" = X (r)dt 2 + 2{f m - / 2 <?i)dMr - (/? - /f)dr 2 - rW (3.16a) 
2M 

X (r) = 1 (3.16b) 

r 

fm - /2ff2 = 1 (3.16c) 

We note that Eq. I|3.16|l specifies the "radial" piece g 11 of the inverse metric only, since the "coordinates" are related to 
the covariant basis elements g^ . Thus, the time part of the gauge is completely free to be chosen, including its direction 
relative to the radial part. This freedom allows for the veritable zoo of forms known for the vacuum metric, which 
can all be obtained from each other via coordinate transformations due to Birkhoff 's theorem. A closer examination 
of Eq. (|3.16() shows that there is only one free function of radius, which is the combination B(r) = f\g% — fagx] once 
this is determined, the other functions display all the possible isometries. 

Particular choices of the fi and gi yield familiar forms for the metric. When fig 2 — f%g\ — 0, the ordinary 
Schwarzschild metric is obtained, while / 2 — = yields the Eddington-Finkelstein metric (advanced and retarded, 
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depending on other sign choices), and ff — /f = 1 gives us the Painleve-Gullstrand line-element, discovered indepen- 
dently by Lasenby, Doran, and Gull p|. For the class of solutions — /| = const., which includes the latter two, see 
^| and the references therein. The Kruskal extension (see e.g. jf|) requires assigning two separate vierbein fields h 
to each point in spacetime, which provides the double coordinate cover. Since this procedure lies somewhat outside 
the scope of this paper, we will not discuss it further here. 

It should be noted that Eq. (|3.1tjf) can be obtained via the usual method of differential forms [Tj], although the 
authors are not aware of any publication of this result. This fact is hardly surprising, since there is not any essentially 
new information contained in this line element. However, we would like to emphasize that the method of geometric 
algebra is completely turn-crank, and, unlike the method of differential forms, the connection coefficients can be 
found directly, as opposed to solving algebraic equations (see 0, for example). In fact, the solution shown here was 
obtained not only by hand calculations, but via a short program written by one of the authors |37| |. 



IV. TRANSLATION OF DIFFERENTIAL FORMS 
A. The Exterior Derivative in Geometric Algebra 

For the sake of comparison with the standard approaches to general relativity (see [H,E3)> it is useful to understand 
how the geometric algebra formalism relates to the differential-geometric formalism. To do this, we will proceed to 
incorporate the language of differential forms into geometric algebra wholesale. 

Before doing anything else, we explicitly identify differential forms with pure multivectors (blades). (Hcstcnes and 
Sobczyk use a slightly different convention, wherein differential forms are scalar-valued 0; see also pjj.) We may 
then expand an arbitrary r-blade in the coordinate basis: 

w = -^ ai a 2 ...a r dx ai A dx a2 A ... A dx ar (4.1) 
r! 

where d is the exterior differential, which we need to understand in the language of geometric algebra. To do so, we 
take the derivative of our arbitrary r-blade, and obtain an r + 1-blade: 



r! \ dx b 

Obviously, dx b is a vector (one-form), while x b is a scalar function. This indicates that d is vector- valued, and acts 
via the exterior product (as expected from its name). Applying the operator a second time yields 

d 2 u = i ( - ^ ai ° 2 ' h ar ) dx c A dx b A dx ai A dx a2 A ... A dx a - = (4.3) 



r! \ dx c dx b 

since the action of partial derivatives commutes, while the wedge product changes sign when two adjacent arguments 
are interchanged. 

Since the coordinate functions x a are scalars, the action of the exterior differential is the same as the gradient is 
the same as the covariant gradient or curl: 

dx a = dx a = Vx a = V A x a = g a . (4.4) 

For generality's sake (and following the precedent set by such texts as we wu ^ consider the covariant gradient. 

In this case, 

w = -.oj aia2 ... ar Vx ai A Vx a2 A ... A Vx a \ (4.5) 
r! 

which is completely equivalent to the previous definition. Applying the covariant curl to this expression gives us 
VAw = if 9u3 T a r ar ) Vx b A Vi" 1 A X7x a2 A ... A Vi° 



r! \ dx b 

+ -r^Ws a I (V A Vx ai ) A Vx a2 A ... A Vx ar - Vi" 1 A (V A Vx a2 ) A ... A Vx a " + . . . 
r! L 

+ (-l) r - 1 Va; Q1 A Vx a2 A ... A (V A Vx ar ) |; (4.6) 
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it is immediately evident that this is equivalent to the exterior derivative if and only if V A Vi" = — which is true 
in the absence of torsion. 

We can extend the development to arbitrary (noncoordinate) frames by transforming the arbitrary multivector with 
a vierbein (frame map) operator: 

7i _1 (w) = -w kl k*...k r fi kl A e k ' 2 A ... A e\ (4.7) 

r! 

where 

e = h k a dx a = h~ l {dx k ). (4.8) 
Applying the exterior derivative to this result, we get 

dh-^uj) = - ( ^-ru> klk2 ... kr h k a \h k a l ■ ■ ■ h k a A dx b A dx a ' A dx a * A ... A dx a % (4.9) 



r 



\dx b 



which follows directly from Eq. (|4.2|) . 



B. The Structure Equations without Torsion 



The orthonormal frame is defined in this way: 

e k ■ e l = r) kl = ~h-\dx k ) -Ti-^dx 1 ) = h^hlsT. (4.10) 
Applying the covariant derivative to this expression yields (in terms of components) 

ei-Ve*=-w^, (4.11) 
which, when protracted, gives us the connection 1-forms: 

V A e k = -e l A (^e J ') = -u> k A e j . (4.12) 

Putting all of this together gives the 1-forms in terms of the connection bivectors: 

u>j h = d a fl(a) ■ (ej A e fe ) . (4.13) 

Although it requires calculating the connection bivectors first, this expression is noteworthy as an explicit expression 
for the connection 1-forms. 

Using the identification of covariant curl with exterior derivative, we get 

de k = V A e k - h^V A Vx^ = -w/ A e? (4.14) 

which we immediately rearrange to yield Cartan's first structure equation: 

de k +ujj k Ae j = 0. (4.15) 

Cartan's second structure equation is found in a similar way: 

doj k + U3i k A Uj l = ®j k (4.16) 

where 

0/ ee R{e 3 A e k ) (4.17) 

is the Riemann curvature tensor 2-form in the orthonormal frame. 

From the preceeding, we see that differential geometry can be folded into the formalism of geometric algebra. In 
addition, Eq. I|4.13|l provides an expression relating the connection 1-forms of differential geometry (for which there 
is no direct expression for calculation that the authors are aware of) to the connection bivectors, which we can find 
explicitly using the techniques of Appendix [C] It remains to be seen whether there is a useful expansion of Eq. (|4.13|) 
in terms of differential forms. 
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C. Example: Connection 1-Forms for the Schwarzschild Case 

On the level of components, the connection 1-forms and bivectors are equivalent, so in a very real sense geometric 
algebra provides a more direct method to computing the connection than is provided by Eq. (|4.15|) . In the absence 
of a more primitive form of Eq. 1|4.13[1 , we use 

u>j k = d a n(a) ■ (e, A e k ) = • (e, Ae')}, 

where the J7 M are the connection bivectors in Eq. (|3.10|1 . to obtain 

<V = Wi° = g[e - g' 2 ei w 2 = u} 2 ° = -5 2 e 2 /r 

^o 3 = ^3° = -92e 3 /r wi 2 = -W2 1 = -gie 2 /r (4.18) 

u>i 3 = —UJ3 1 = —gies/r W2 3 = — W3 2 = — cotOe^/r 

These results reduce to the canonical expressions (see e.g. 0) used in obtaining the Schwarzschild solution when 
g-2 = 0, and agree with the expressions obtained using differential geometry directly [14j. 

V. NULL VECTORS AND TETRADS 

Null vectors, especially null geodesies, often indicate symmetries in particular spacetime geometries. In standard 
expositions |l8l Il3 . |2(1 l2l| , a complex structure is often introduced to allow for four null vectors in a spacetime, and 
tensorial objects acquire complex character as well. This is known as the null tetrad formalism, and is intimately 
related to spinors. 

Although we will not discuss spinors in this paper, null vectors and tetrads are useful to our general development. 
An external complex structure is not needed in geometric algebra, since we have an object ready-made to fill the 
role of i: the spacetime pseudoscalar. We will "complexify" spacetime in the next section, then define null tetrads as 
mixed-grade multivectors in Section IV Bl which geometrically are Penrose's flags and flagpoles. In Section lvTl we will 
find that the "complex Weyl tensor" is identical to the Weyl tensor we have already defined in Section [H] and discuss 
the Petrov classification. In this way, geometric algebra encompasses the complex structure often used in general 
relativity, elucidating the underlying geometric structures, without introducing any fundamentally new objects. 

A. Paravectors and "Complex" 4-Space 

As discussed in some detail in Appendix IA 31 the space of scalars, pseudoscalars, and bivectors form an algebra 
under the geometric product. This is called the even subalgebra C£f 3 , and it is well-known that it is isomorphic to 
the three-dimensional geometric algebra CI3. We will find it useful in dealing with the Lorentz group @], and in 
constructing a (3 + l)-dimensional "complex" space [22,|23|]. The latter discussion will allow us to construct null 
tetrads, which in turn leads us to the Petrov classification of the Weyl tensor, which we will examine in Section IVll 
In a subsequent paper, we will examine the topic of 2-spinors in this space; see also 24] for a preliminary sketch of 
some of these ideas. 

A useful way to obtain the even subalgebra requires selecting a particular timelike direction t in the full algebra 
C£i : 3, where i = 1 for simplicity This allows us to define a paravector a for any vector a in spacetime: 

a — > a = at = clq + A (5-1) 

where ao = a ■ t and A = a A t define the scalar and timelike bivector parts of the paravector. As discussed in 
Appendix I A 31 the space formed from scalars, pseudoscalars, and bivectors is closed under the geometric product, so 
the product of paravectors with paravectors will always lie in this subalgebra. 
If we define the operation of Hermitian conjugation for any multivector Q as 

= tQt, (5.2) 

we see that the paravector we defined in Eq. H5.1JI is Hermitian: a} = a. In fact, the selection of the timelike direction 
allows us to split the even subalgebra into two pieces: Hermitian, or "real", and anti-Hermitian, or "imaginary". In 
terms of a generic multivector tp S Cl\ 3 



(5.3) 
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this split can be written as 
where 

(±)z = (5 - 5) 

In the remainder of this paper, we will often talk about "real" and "imaginary" bivectors (meaning timelike and 
spacelike), and use the term "complex number" to mean the sum of a scalar and a pseudoscalar. 

For simplicity, let us define a (real) basis for C£f 3 that contains the elements {l,e = <ti,u = <T2,v = 0-3}, where 
e 2 = u 2 = v 2 = 1 are real bivectors. We also require that the bivector basis be right-handed: eu = — ue = iv and 
cyclic permutations, so that 

a- l cr j = tn ■ <Tj + tTi x a-j, (5.6) 

which is to say that the outer product of real bivectors vanishes. In this basis, a "complex paravector" is written 

V> = "00 + i> e e + ip u u + ipvV, (5-7) 

where the ip a are "complex" scalars. A complex paravector can therefore be seen as a four-dimensional complex 
vector, whose basis elements are the scalar and three real bivector basis elements, and whose components are the four 
complex (or eight real) scalars. (Note that a complex paravector is isomorphic to a Dirac spinor in the language of 
Hestenes 251; see also 0, Chapter 10.) 

Because paravectors are mixed-grade multivectors, the usual inner and outer products are often of limited usefulness; 
instead, we define (echoing Baylis |23) "scalar" and "bivector" products: 

(AB) S = (BA) S = \ (AB + BA) (5.8a) 
(AB) B = -(BA) B = i(AB-BA) (5.8b) 

where the result of Eq. H5.8afl is a complex number. (Recall that A denotes the reverse of a multivector.) In the case 
of a "real" paravector, the scalar product returns the Minkowski metric: 

(ab) s = aobo — a • b = a ■ b. (5-9) 

Thus, paravectors can be used as a replacement for spacetime 4-vectors; for a formulation of relativistic electromag- 
netism entirely in the even subalgebra C£f 3 ~ C£ 3 , see [22^ . 



B. Null Tetrads 



In physical applications, null vectors often play an important role: null geodesies identify paths for light, and can 
indicate symmetries in a system. They are often used in obtaining the Petrov classification of the Weyl tensor, and 
are used in the Newman-Penrose method of finding exact solutions to Einstein's equations 19], which we will not 
discuss in this paper. A particularly useful way to work with null vectors is through complex null tetrads, which 
draws on the previous section. 

If we let TV 2 = be an arbitrary null vector, the paravector 

N = Nt = N (l + r) (5.10) 

places all the relevant information in two quantities: the "extent" No = N • t, and the spatial direction r, which is 
defined to be of unit length. It is easily seen that ( NN ) S = 0, following from Eq. I|5.9|l . For the sake of discussion, 
we arbitrarily align r with e of the paravector basis from the previous section, and let 

7V = 27V [ where Z = -(l + e). (5.11) 
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The paravector I is called a flagpole, which is intimately related to spinors (see [24J for more information). It has a 
few other interesting properties, some of which we summarize here: 

i 2 = I idempotence 1 + 1=1 complementarity (5 12) 

U = nullity t} =1 Hcrmicity ^ ' ' 

These properties indicate that I and I are linearly independent, so we define a new basis for paravectors around them: 

l = ^(1 + e) 

n = 1=1(1 ~e) 

1 . . 1 . 
221 = -(1 + e)u = -(u + jv) 

mt = I u (l+e) = i(u-iv). (5.13) 

We call this set a rawH tetrad, since its elements obey the following relation: 

11 = nn = mm = m [ rh? = (5-14) 

and 

( ln )s = - ( mmt )s = \ ( 5 - 15 ) 

with all other scalar products vanishing. The first two tetrad elements in Eq. (|5.13|) are flagpoles, while the second 
two are their corres pon ding flags. (Since we do not deal with spinors in this paper, the construction of tetrads is a 
trifle awkward. See \l \ for a better technique in geometric algebra, and [2lJ for the canonical discussion.) 

We will occasionally find it useful to let E_ = I, E_i = n, E_ 2 = ZZii an( i ^3 = 221^ i m this notation, we see that the 
basis transforms as 

= RE a R\ (5.16) 

where R is a rotor, as before. We examine the transformation that preserves a particular null direction in Appendix iDl 
The null tetrad "metric" , which is invariant under Lorentz transformations, is 



2(M a E b )s = Vab (5.17) 



where LU 



/o 1 
10 
-1 

Vo -1 



(5.18) 



A general multivector is then written as 



A = A a E a , (5.19) 



where the A a are complex scalars — demonstrating that the null tetrads form a basis for C^ 3 over the complex 
numbers. 

The bivector basis of |18j (with different indexing, however) can be reproduced by taking the nonzero vector products 
of the tetrads: 

Z = -V2(lm) E = -=(1 + e)u 
v2 

Zi = (In) B — \ mrn} ) g = e 

Z 9 = -y^(nmt) = -L(l- e )u (5.20) 
v2 
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These are linearly-independent, and span a three-dimensional space over the complexes. We can write a "metric" for 
them: 



(Z a Z b ) s = z ab -> [z ab ] =0-10 (5.21) 




which is obviously real-valued, despite the "complex" character of the Z a themselves. Thus, when we write a bivector 

A = A a Z a (5.22) 

the complex components A a are easily extracted using the scalar product. And, unlike in the standard formalism, 
the bivector basis is easily interpreted in geometric algebra: Z and Z 2 are oppositely-directed null flags determined 
by the timelike direction Zi = e. We will see further examples of the new geometric interpretations that geometric 
algebra provides in the next section. 



VI. CLASSIFICATION OF THE WEYL TENSOR 



This section is concerned not with finding solutions to the field equations, but rather extracting geometrical in- 
formation about the Weyl conformal tensor. In the case of vacuum solutions, the Weyl tensor is the gravitational 
field, so being able to discuss the form of the solution without first solving the equations can help us understand the 
underlying physical configuration. The classification of the Weyl tensor due to Petrov, which has been translated 
into the language of geometric algebra by Sobczyk 0, l23l l2t| and Hestenes and Sobczyk Q , provides one such way 
of describing the fields, independent of a coordinate system. We examine an algebraic technique using the eigenvalue 
problem in Section fVI Al then we apply the null tetrad formalism in Section fVI Bl to provide a geometric interpretation 
of certain of these classes. This last discussion will result in illuminating new canonical forms for the gravitational 
field. 

The Weyl tensor is symmetric 

A • C(B) = C(A) • B (6.1) 

and self-dual 

C(iB)=iC(B), (6.2) 

where A and B are bivectors. Combining these properties yields the important result £| 

(AC(B)) S = (BC(A)) S . (6.3) 

We draw on an electromagnetic analogue by writing |l8j | 

C(B) = £(B) + iS(B) (6.4) 

where 6(B) and 6(B) are both real bivector-valued symmetric operators. Because of the self-duality of the full Weyl 
tensor, we see that 

E(iB) = —6(B) and 6(iB) = E(B). (6.5) 
These results will be used extensively in the following sections. 



A. The Petrov Classification: The Eigenvalue Problem 



The Petrov type is most easily determined by solving the eigenvalue problem: 

C(B) — AB = (6.6) 

where A is a complex scalar. The Weyl tensor can be viewed as a 3 x 3 complex symmetric matrix acting on the 
3-dimensional complex space of bivectors. Since the operator C — XI is zero, it is singular, and we are left with the 
usual determinant 



[C — Al| =0. 



(6.7) 
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If we define a timelike (real) bivector basis {cfk} (for k = 1,2, 3), we can write the determinant as 

\C-\1\ = -ffio-aira <[(C(o-i) - Ao-i) x (C(<x 2 )-A<t 2 )](C(<x 3 )-A<t 3 )) s = (6.8) 

which is a third-order algebraic equation in A. Although this expression looks hideously complicated, we will shortly 
calculate it for the Schwarzschild case, and show that it is actually remarkably simple to work out. (See also (2^ for 
a worked example in C£ 3 .) Because the Weyl tensor is tractionless 

d a C{aAb) = Vb (6.9) 

the eigenvalues satisfy the condition 

Ai + A 2 + A 3 = 0, (6.10) 

which means only two eigenvalues are independent. 

Stated simply, the eigenbivectors determine the general Petrov type (labeled I, II, and III), while specific eigenvalue 
restrictions yield special cases (types D and N). Petrov type I occurs when the eigenbivectors span the entire 3- 
dimensional space; the bivectors for types II and III span a two- and one-dimensional subspace, respectively, and both 
have a null eigvenbivector: 

Bl = (types II, N, and III) 

The last two types are special cases: D corresponds to type I when two eigenvalues arc equal, while type N corresponds 
to type II with all three eigenvalues vanishing. Geometric interpretations are left to Section fVl Bl 
Hestenes and Sobczyk 3] have shown that the Weyl tensor can be written as 

C(X) = Q jk ajXcr k (6.11) 

where the symmetric matrix Qi — Q 3 is comprised of complex scalars. However, with the eigenbivectors in hand 
for types I and II (which include types D and N), we can rewrite this as 

C(X) = £j jB fc XB fc + ixB£j// fc (types I, II, D, and N) (6.12) 

where fi% = Ai + A 2 /2 and /i 2 = A 2 + Ai/2 are linear combinations of the eigenvalues. (For type I, the third eigenvalue 
is found from Eq. I|6.1U|I . and the third eigenvector is not needed for the calculation.) Type D corresponds either to 
Mi = M2 01 M2 = 0, while type N corresponds to fa = when is the null eigenbivector. Type III has the canonical 
form 

C(X) = (BXF + FXB) ft (type III) (6.13) 
where /i is a complex scalar, and the bivector F is orthogonal to the eigenbivector: F • B = 0. 



B. Geometric Classification Using Tetrad Techniques 

The forms for the Weyl tensor obtained in the previous section are useful, but we can go a step further. Motivated 
by the Newman-Penrose formulation, we use the tetrads and the bivector basis constructed from them to understand 
what the types mean geometrically. Although Szekeres obtained the same qualitative results psj . here we provide 
new canonical forms for the Weyl tensor. 

We begin by using Eq. I|6.3|) to define the complex scalar-valued function 

(A, B) = - (AC(B)) S = A ■ C(B) + A A C(B). (6.14) 
When we use the bivector basis defined in Eq. I|5.20[l . we get 



* Qb = -(Z Q C(Z b )) s = * b , 



(6.15) 
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*0 


= *00 




*1 


= #01 


= * 




= * n 


-#< 


#3 


= lj 12 


= * 


*4 


= #22 





which can be regarded as the components of a 3 x 3 symmetric, complex matrix. We reduce the amount of redundancy 
by labeling the components |19| 



2 = * u - *2o = - (Z 1 C(Z 1 )) S + <Z C(Z 2 )) S 
!1 = -(ZiC(Z a )) s 

(Z 2 C(Z 2 )) S . (6.16) 

The sixth component, which has been folded into # 2 , 

#02 = #20 = -(Z C(Z 2 )) S (6.17) 

is generally not talked about in algebraic discussions of the Weyl tensor; this is because it is not an independent 
quantity, due to the trace-free property jig. In certain cases, it is useful to fix one null direction and rotate the 
others, thereby obtaining new values for some of the # a ; Appendix IDl is devoted to that topic. 

To perform certain calculations, we will find it useful to write the complex bivector basis {Z a } in terms of the real 
bivector basis {e, u, v}: 

e = Zi u=4=(Zo + Z 2 ) v = —=(Z -Z 2 ), (6.18) 
V2 i\ 2 

and we split the tensor into its electric and magnetic parts and consider the quantities 

* ob = -(Z Q C(Z b )) s = 2- 1 / 2 (E u (Z b )-B v (Z b )+i{E v (Z b ) + B u (Z b )}) (6.19a) 
*16 = -(Z 1 C(Z b )) s = E e (Z b )+iB e (Z b ) (6.19b) 
# 26 = -(Z 2 C(Z b )) s = 2- 1 / 2 {E u (Z b ) + B v (Z b )-i{E v {Z b )~B u {Z b )}) (6.19c) 

along with the symmetry conditions 

* 01 - * 10 ^ 1 V2B e {z ) = ^( Zl ) + b.czo (6 ' 20a) 

, Tf _. T . / E U (Z 2 ) - B V (Z 2 ) = E U (Z ) + B V (Z ) (Ronu\ 

* 02 -* 20 ^ 1 E V (Z 2 ) + B U (Z 2 ) = -E V (Z ) + B U (Z ) ( 6 - 20b ) 

0, J v / 2 J B e (Z 2 ) = £ u (Z 1 )+S K (Z 1 ) 

* 12 - * 21 ^ { V25 e (Z 2 ) = -^( Zl ) + B„( Zl ) (6 ' 20C) 

In the preceeding, i? e (Zb) = e • E(Z b ), and so forth. 
We will use these decompositions to obtain the Petrov types for the various # a in the following subsections. 

1. Transverse Waves: \&o 7^ or ^4 7^ 
In the case where $0^0 but the rest vanish, we immediately get (from Eq. Ij6.19|) ) 



while applying Eq. 1|6.20J) yields 



and 



E e (Z b ) = B e (Z b ) = 0, 
C(Zi) = C(Z 2 ) =0 



B V (Z ) — —E U (Z ) 
E V (Z ) = B U (Z ), 



so the electric and magnetic pieces are perpendicular to each other. Putting all of this together, we see that 

*o = V2{E u (Z )-iB u (Z )) 
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and 

E(Z ) • 6(Z ) = E U {Z )B U (Z ) + E V (Z )B V (Z ) = E U {Z )B U {Z ) - B U (Z Q )E U (Z ) = 0. 
We can then write the nonvanishing piece of the Weyl tensor as 

C(Z ) = E U Z 2 = -B V Z 2 , * ^0, 
or, utilizing Eq. 15.21fl . for an arbitrary bivector X = X a Z a 

C(X) = -£ U Z 2 (Z 2 X) S = ^Z 2 XZ 2 . (6.21) 

Similarly, when ^4 7^ but the rest vanishing we get 

* Q#4 = -> C(X) = ^Z XZ , C(Z ) = C(Zx) - 0. (6.22) 

Equations l|6.21|) and (|6.22() already are in the canonical form for Petrov type N; confirming this classification is a 
simple matter, using Eq. (|6.8|l in the form 

\C-X1\ = -euv([(C(e) - Ae) x (C(u) - Au)] (C(v) - Av)) s = 0, (6.23) 

and using the canonical forms to evaluate the expression. 

The Weyl tensor represented in both cases here is proportional to a null bivector. Its electric and magnetic pieces 
are perpendicular to each other, and both are perpendicular to the direction =pe. This means that the Weyl tensor 
represents a transverse plane wave traveling in the (=Fe)-direction, which agrees with the analysis by Szekeres |2S| 
(albeit by a different route); see for the electromagnetic analogue. 

2. Coulombic Case: <3> 2 / 

The final case which is directly related to the electromagnetic treatment is when ^ = 0. Here, we find immedi- 
ately that 

C(Zi) = (£ e (Zi) + iB e {Z x )) Zj = F 1 Z 1) (6.24) 

while 

C(Z ) = (E u (Z )+iB u {Z ))Z = F 2 Z 
C(Z 2 ) = (E u (Z )+iB u (Z ))Z 2 =F 2 Z 2 , 

where we have exploited the symmetry conditions to get 

E e (Z ) = E e (Z 2 ) =B e (Z ) = B e (Z 2 ) = 
B B (Zi)=^(Zi)=fl„(Zi) = B (Zi) = 

and 

E U {Z 2 ) = E u {Zq) B u (Z 2 ) = B U (Z ). 

The complex scalars F a are obviously the eigenvalues for the system, which we can confirm by putting all of this 
into Eq. (HT23J): 

Ai = E e {Z x ) + iSe(Zi) =F U A 2 = A 3 = E U {Z ) + iB u {Z ) = F 2 . 
Enforcing Eq. I|6.10(l further restricts the eigenvalues: 

3 1 1 

VA, =0 -» E a (Z ) = -~E e {Z x ) and B U (Z ) = --B e (Z x ). 



k=l 
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This is Petrov type D, with eigenvectors Z& (or equivalently {e,u, v}). This means that all incoming information is 
balanced by outgoing, if one regards e as a "radial" direction — hence, type D tensors are "Coulombic" |20j . 
We obtain the canonical form by noting that 

C(X) = -F 1 Z 1 (XZi) s - F 2 (Z (XZ 2 ) S + Z 2 (XZ ) S ) = ±Fi |z 1 XZ x - ~ (Z XZ 2 + Z 2 XZ ) 

Since Zg = Z 2 = and Zi anticommutes with both Zq and Z 2 , we discover the identity 

i (Z XZ 2 + Z 2 XZ ) = -X aX -S aX -L x (Z Z 2 + Z 2 Z ) = (XZ^ Z 1; 

so that 

C(x) - i^i {ZiXZ! - (xz 1 ) s zj = ~Fj jziXZi + ~x j , (6.25) 

which is the canonical form for Petrov type D. 



3. Longitudinal Waves: $1 /fl or ^3 7^ 

For *i ^ 0, we find that 

*i = E e {Z ) + iB e {Z ) 

using the same techniques as before, while the rest are zero; the Weyl tensor is then 

C(Z ) = (E e (Z )+iB e (Z ))Z 1 = FZ 1 
C(Z0 = (£ e (Z ) + z5 e (Z ))Z 2 = FZ 2 

C(Z 2 ) = 0. (6.26) 

C(Zp) is solely in the Zi-direction, while C{Z\) is a transverse piece — so this object has a longitudinal component 
[28j. The characteristic equation gives us three zero eigenvalues, and since Z 2 is the only eigenbivector, we have 
Petrov type III. The canonical form is found immediately to agree with that analysis: 

C(X) = —F ((XZ 2 ) S Z x + (XZi) s Z 2 ) = I (Z 2 XZ! + Z!XZ 2 ) . (6.27) 
Repeating the process for ^3 7^ is trivial. 



C. Example: The Petrov Type for the Schwarzschild Case 



For simplicity, we will use the orthonormal frame from Section 1111 ij| to define a natural bivector basis for the 
Schwarzschild solution: 

where k runs from 1 to 3, and cr 1 cr 2 <T3 = i. In this frame, the curvature bivectors are 

/?( CTl ) = C{ax) = -ix"(r)o-i , /?(<x 2 ) = C(<r 2 ) = ^(rjffj , /?(<x 3 ) = C(<r 3 ) = ^x'(r> 3 

where we have used the vacuum equations to set the Riemann tensor equal to the Weyl tensor; recall that x( r ) = 
1 — 2M/r. The eigenvalues and eigenvectors of C can be read off immediately; however, for the sake of example we 
will evaluate Eq. I|6.8|l to obtain them. 

The characteristic equation requires evaluation of the quantity 

(C(<ri) - Ao-x) x (C(<7 2 ) - A<r 2 ) = (-\x"(r) - \ \ (^x'(r) — \ \ <T\ X ct 2 ; 
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noting that cr\ x <r 2 = icr^, we get 

\C-X1\ = -i (-^"(r) - a) ^ X '(r) - \) (i<r 3 <r 3 ) s = (-^x"(r) - a) (^x'(r) - ^ ' 
Solving this gives the eigenvalues 

x 1 »M 2M x x 1 v x M 

two are repeated, so this is automatically either type D or (if \ is constant) type 0, which is the trivial Minkowski 
case. The eigenbivectors are 

Bfe = cr k , 

which spans the full three-dimensional complex space. 

From Eq. (|(j.25[) . we see that the canonical form of the Weyl tensor is [BJ 

C(X) = -|^( ( x 1 X<x 1 + V) (6.28) 

since fix — Ai + 1/2A 2 = — (3M)/(2r 3 ) and ^ 2 = A 2 + l/2Ai = 0, and <j\ = 1. The eigenvalues are real, so there is 
no "magnetic" piece: the nonrotating black hole solution is purely "electric". As we will see in Section lYlII the Kerr 
solution, which describes a rotating black hole, does contain a magnetic part, and in fact a trick exists to obtain it 
from the Schwarzschild solution by replacing the radial coordinate r with r — ia cos 9 0, |29| . 

VII. EXAMPLE: THE KERR SOLUTION 

The Schwarzschild metric is among the simplest solutions to Einstein's equations, but physical black holes will 
generally have nonvanishing angular momentum. A more realistic — and more complicated — solution is the Kerr 
metric. No truly simple techniques exist to obtain this metric, but a straightfoward method has been developed by 
Schiffer et al. 30] , and adapted to the language of geometric algebra by Doran and colleagues |2^, , along with 
Snygg . We will follow their path, showing the power of geometric algebra by performing coordinate- independent 
calculations, while providing some physical insight into the meaning of the ansatz used. 

A. The Kerr-Schild Ansatz and Curvature 

The Kerr-Schild ansatz, which is used by [H IH III El H2 , 

is defined in vierbein language as 

h~ 1 {a)=a + MVa-ll, (7.1) 

where I 2 = is an unspecified null vector, M is identified with the black hole mass, and V is a scalar function. From 
this ansatz, we can define a new frame and its reciprocal from an arbitrary orthonormal frame 

5 M = iT 1 ^) =e ll +MVe tl -ll 

= Tj^) = e M - MVe^ ■ 11. (7.2) 

From Eq. I|7.2|) . we obtain the usual Kerr-Schild metric 

9V =9^gu = V»v + 2MVe.fi ■ le v ■ I. (7.3) 

Two other features of the KS vierbein will prove to be important for us: 

h(l) =hr 1 (l) =1 and /? _1 (a)i = h(a)l = al, (7.4) 

which means that I has the same components in both frames — 

I = g^g^ ■ I = e% • I 
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— while 

= a ■ V(Z ■ I) = 21 ■ (a ■ Vl) . (7.5) 

This ansatz is also useful for simplifying the calculation of connection bivectors: using Eq. (SB) in Appendix O we 
see that the third term vanishes, while the first two terms can be rewritten as 

ft (6) = -M (V A (VI ■ bl) - Vg x A 11 ■ (v A &)) • (7.6) 

From this and Eq. I|7.5|l . we see immediately that ft(Z) = 0, while if b = g a or b — e a , the second term vanishes: 

ft,, = Sl( 9fi ) = n(e M ) = -MV A (VIJ), (7.7) 

using the shorthand 1^ = • I = g„ ■ I. If we wish to write everything out in terms of the orthonormal frame, we get 
another interesting expression: 

ft (6) = M {MV 2 l ■ bl A v + Ve x A 11 ■ (v A b) - e A A V A (VI ■ W)} (7-8) 

where 

v = l-Vl = l-Vl (7.9) 

is a vector that vanishes when I is geodesic; note that I ■ v = via Eq. H7.5(l . 

Due to the simple form of the connection bivector, we can write down the curvature tensor involving I immediately: 

R(g„M) = -l-Vn fi (7.10) 

Following |3l| and [32l | , we form the quantity 

l-R(l) = l-(g^R( 9ll M)) (7.11) 

where 

g» ■ R(g^ A I) = -g» ■ [l ■ Vft M } = -I ■ V (g» ■ ft M ) + (l ■ V g ^) • ft M (7.12) 

is the Ricci tensor R(l). Noting that 

- M" V ■ ft /j = + ZV • (VI) = Vv + I (l ■ VV + VV ■ l\ (7.13a) 

-M~H ■ Vg^ = I ■ V (VW) = VvF + I (in ■ VV + Ve^ ■ v^j , (7.13b) 

and that ^ M ft M = ft(Z) = 0, we find 

A/" 1 /?(Z) = 2vl ■ W + VI ■ Vv + vVV l + U Ve-Vl- Sl(v), (7.14) 

where 

e = v - (vi) = i- vv + vv -i 

is the expansion parameter for the vector VI. Taking the inner product of Eq. 17.14JI with Z eliminates all but one 
term: 

Z • R(l) = MVl ■ (l ■ Vv] = —MVv 2 , (7.15) 



where the last expression follows from the vanishing of v ■ I. When v vanishes, I is a null geodesic, and the Ricci 
tensor satisfies Einstein's equations. (Note that any vector v cx I also makes Eq. I|7.15|l vanish, but we can transform 
the proportionality to zero simply by redefining the time coordinate 1 0J . ) From here on, we let 

v = 0, 
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even when considering Einstein's equations in the presence of nonvanishing energy-momentum. 
Plugging the geodesic condition into Eq. (|7.14|) yields an eigenvalue equation: 

r(i) = imi ■ ve = im{(i ■ v) 2 y + {i ■ w)v • i + vi ■ v(v (7.16) 

so that the vacuum condition is I ■ V0 = 0. The connection bivector becomes 

ft M = -Me x AV A (y;g. (7.17) 

Using this knowledge, we determine that I is a principal null direction by forming the quantity 
I ■ C(a A I) = an ■ A I) - 1 ( ff „ A /?(/) + Z A + iiZ ff|t A «J 

where 

l-R(g / ,Al)=Mll„ (Vv)V 

i • A R(l) + I A = I ■ g„R(l) - 11 ■ R(g^) = 

(using the symmetry of the Ricci tensor I ■ R(g^) = g^ ■ R{1))- What remains is 

I ■ C{a M) = Hll ■ a (7.18) 

which is an eigenvalue equation. Obviously, in a vacuum the Ricci scalar vanishes, but this does not change the 
essence of the Weyl tensor. Along the same lines is the expression 

I A C{a A I) = I A R(a A I) = Ml ■ a (l ■ w) I A V A I. (7.19) 

From this development, we see that I ■ (a A I) = 11 ■ a (and hence I) is a principal null direction, which we will examine 
in more detail in the subsequent section. We will confine ourselves to the vacuum case from here on. 

B. Kerr Geometry 

To recap, we have found that I of the Kerr-Schild ansatz is a principal null direction. From it, we define the bivector 

e = I A e , (7.20) 

and write 

Z = (l + e)e (7.21) 
so that lo = I ■ g Q = 1 and e 2 = 1. From this null vector we define the tetrad paravector 

[=^e = i(l + e) (7.22) 

and select an arbitrary unit vector that satisfies the properties 

u 2 = — 1 , u ■ I = , u ■ eo = 

to generate the rest of the tetrad. More important to us for the moment is the bivector basis from Eq. I|5.20|l . which 
we write as 

V2Z = uAl Zi=e V2Z 2 = u A (e ie ) . (7.23) 
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From these, we can form the <J/ a as in Section IYl Bl The first two are identically zero, following from Eq. (|7.18(l and 
Eq. fTHty : 

*o = ^ {(« A I) ■ C(u A I) + (ii A Z) A C(« A 0} = ^ {« • (I ■ C(« A Z)) + u A (Z A C(« A I))} = 

*i = ~ {(e A I) ■ C(u A I) + (e A Z) A C(« A I)} = i {e • (Z • C(« A Z)) + e A (Z A C(« A Z))} = 0, 

since u ■ I = by definition. This demonstrates that any Kerr-Schild solution is algebraically special. The actual 
Petrov type, however, depends on other choices, which we will make now. 

To obtain a black-hole solution, we require that V and e be functions of the spatial coordinates only. With this 
criterion, we will reduce the field equations to expressions in the three dimensional subalgebra C£i 3 . Since VoA = 0, 

we define the Cartesian 3-gradient V= e fc Vfc and its bivector equivalent 

V = e A V= e V . (7.24) 
Thus, the four-dimensional gradient can be reduced to 

V = g k V k =V -M Vie ■ V, (7.25) 

from which we obtain I ■ V = e • V. 

Following |3Cj . we now obtain the field equations; since M is an arbitrary constant, each coefficient of a power of 
M must vanish separately. The Ricci tensor is 



Ru = g^- R{g^ A g v ) = g^ ■ (v M ft„ - V„ft M + ft M x ; 



each $1 contributes a factor of M, while g 11 = e M — MVIH contains a zero-order and first-order piece. Therefore, we 
have the following expressions: 



0{M) : e» ■ (V M a - V v f2 M J -V 0, - V, (e^ • H M ) (7.26a) 

0{M 2 ) : -MVIH ■ - VMn) + e M ■ (fl M x ft u ) (7.26b) 

C(M 3 ) : -MVIH ■ (SI,, x U v ) = -MVl ■ (Sl(l) x = 0. (7.26c) 

As the reader can easily tell, we have reduced the complexity of the problem significantly: the first-order equation, 
which we will treat momentarily, is straightforward to find, while the third-order expression is an identity. The 
remaining equation is satisfied by the solution to the first-order equation; we will not prove this here, since it does 
not involve anything new. 

0(M) Field Equations 

After some work we find that 

m- 1 v n v =v 2 (vu„) + v (z„e) , 

where 8 = V • (VI) — e ■ 'VV + W ■ e as before, and we have used the vanishing of I ■ V0 = e ■ V9 from Eq. I|7.16|l . 
We already know that 

■ o M = -g» ■ = Mie, 

so that the first-order field equations are 

V 2 (^ZZ,)+V(Z,e) + V,(Z9) = 0. (7.27) 

When v = 0, Eq. lf7"27f) reduces to 

V 2 (VI) + V 6 = V 2 {V (e + e • e ) } + e • V9 = 0; (7.28) 
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multiplying by <=o splits the equation into scalar and bivector pieces: 

V 2 V = Q and V 2 (Ve) + V6 = 0. (7.29) 
For the case when v — k (for k — 1,2, 3), we have 

V 2 (Vll k ) + V (l k B) + V fe (19) = | V 2 (Veoh) + e V fc e| + e> V 2 (FZ,-Z fc ) + V (i fc e) + V fc (e9) • e ; (7.30) 

the term in braces vanishes due to Eq. (|7.29() . Expanding the Laplacian yields 



V 1 (Vljh) = l,V z (Vl k ) + l k V z (Vlj) + 2V \V IA ■ \V l k j = ljV k Q + hVjQ + 2V \V IA ■ \ V l k j (7.31) 
so that the field equations are 

2Ve 3 (v ■ (v l^j + 6 jv l k + (V fc e) ■ e | = (7.32) 

or, in bivector language, 

2V {V (eo A e k ) ■ e} ■ Ve + 9 [V {(e A e k ) ■ e} + (e A e k ) • Ve] = 0. (7.33) 

Our task is now to find a reasonable expression for the derivatives of e; once again, we follow |30j . 
For some arbitrary vector a we define the bivector a = ^ A a, and from it create the operator 

M(a) = a • Ve, (7.34) 

whose adjoint is 

M(a) = V (a • e) = V (a • e) - d h ((b • Va) • e) ; (7.35) 

we have used the accents to show that a is not differentiated. From these definitions, we see that the field equations 
can be rewritten as 

^MM(a) = M(a) + M(a), (7.36) 

and we note that M(e) = M(e) = 0. Due to the latter property, we make the following ansatz: 

M(a) = a U(a) (7.37) 

where U_ is a unitary operator, which obeys U_{e) = U(e) = e. A unitary operator that preserves a certain spatial 
direction can be written as a rotation: 

U(a) = e^^ae-"^/ 2 (7.38) 
which means that, after expanding the exponentials, we find 

a • Ve = M(a) = a(a x e)e + (3i(a x e), (7.39) 

where 

a = (1 — cos ^i) and /3 = sin-0. (7.40) 

(Recall that the symbol x represents the commutator product; if we were using the vector cross product, the pseu- 

doscalar element i would not be present.) Removing the arbitrary bivector a gives us the divergence and curl: 

Ve = <9 a a • Ve = 2a + 2(3ie 

-> V • e = 2a and V x e = 2if3e. (7-41) 
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Following [2|J again, we form the Laplacian of e via Eq. I|7.39[l and by direct calculation (in order of mention): 

V 2 e = V • Ve = Va - (e • Va) e - 2 (a 2 + /3 2 ) e + ie x V/3 (7.42a) 
V 2 e = (V(Ve)> 2 = 2 (Va -iVx (/3e)) = 2 (Va + ie x V/3 - 2/3 2 e) (7.42b) 

where we have used the identity V x (/3e) = /3V x e — e x V/3. If we further note that 

V • {/3b) = = /3V ■ e + e ■ V/3 = 2a/3 + e ■ V/? 

and combine the Laplacian expressions, we obtain 



Va = (/3 2 - a 2 )e + ie x V/3 V/3 = -2a/3e - ie x Va 
e • Va = /3 2 - a 2 e • V/3 = -2a/3. 



(7.43) 



From this, a little more work shows that 



and 



If we let 7 = a + i/3, we can collect the directional derivatives into a single equation: 

e • V 7 = - 7 2 , (7.44) 

while 

V7 = -7 2 e + e x V7. (7.45) 

V 2 7 = 0, (7.46) 

(V7- 1 ) 2 = 1; (7.47) 

thus, 7 is harmonic, and its inverse satisfies the eikonal equation. 

Since both V and a are harmonic, it is reasonable to ask whether they are equal; it turns out this guess is consistent, 
by evaluating Eq. I|7.29[) . Noting that V — a yields 

9 = e- Va + aV-e = a 2 +/3 2 = 77 f , (7.48) 

we work out the Laplacian using Eq. (|7.42(l and Eq. (|7.43|) to give us 

V 2 (ae) = aV 2 e + 2(Va) • Ve = 2(aVa + /3V/3) = V(a 2 + /3 2 ). (7.49) 

With all of this in hand, we find we can write e in terms of 7 : 

i + (v 7 - 1 )-(v( 7 - 1 )t) v ' 

The Schwarzschild solution is found by letting 7 = V = 1/r, so that j3 = 0. The Kerr solution is a little more 
complicated: 



7~ 2 = x 2 + y z + (z - ia,y 
V = <7> 



p 4 + a 2 z 2 

1 z 

((px + ay)crx + {py - ax)cr 2 ) + -cr 3 (7-51) 



,2 1 n 2 



where p is implicitly defined by 

4 2/ 2 2\22 n 2_2,2,2 in rr,\ 

p— p(r— a) — az=Q; r = x +y + z . (7.52) 

This value for 7 is picked to ensure that the system is axisymmetric (about the z-axis); the parameter a is related to 
the angular momentum of the black hole. 
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VIII. FURTHER DIRECTIONS 

We have by no means exhausted the task of translating solution techniques into the geometric algebra language; 
indeed, we believe that further insights into the geometrical structure of general relativity can be obtained through 
this process, and perhaps more techniques can be discovered. In addition, we point to a paper by Sobczyk 33], 
which discusses the embedding of the spacetime of general relativity in a pseudo-Euclidean space of higher dimension, 
utilizing geometric algebra throughout. 

One area of interest is the extension of these techniques to the space of Lorentz spinors (also known as 2-spinors 
The work in Section l5.13l is a hint toward this, as are the papers by Jones and Baylis J0, and Lasenby, Doran, 
and Gull [34| ■ Spinor techniques lead to an adaptation of the Newman-Penrose formalism |l9j|. which has been useful 
for obtaining many exact solutions. 

The viewpoint in the first section of this paper, in which an orthonormal frame of vectors is considered to be 
fundamental (as opposed to taking a coordinate metric as fundamental) gives a somewhat different ontological status 
of gravitation, compared with the usual interpretation of general relativity. Specifically, this construction can be 
interpreted as a flat spacetime theory, in which gravitation is a modification of the geometric relationships of an 
underlying flat spacetime; mathematically, this is locally equivalent to general relativity, which is commonly viewed as a 
theory of intrinsically curved spacetime. We also emphasize that these two viewpoints are physically indistinguishable, 
as geodesies are the same in both. 

If gravitation is viewed as occurring in a Minkowski space arena, it becomes conceptually closer to the other 
fundamental forces. Indeed, the same covariant derivative can also be obtained in a flat background spacetime by 
enforcing local gauge invariance under the gauge group SL(2,C), the proper Lorentz group [5|,l3- Many fundamental 
questions related to these gauge formulations remain unanswered, but the existence of these theories stands as an 
interesting contrast to the usual approaches to gravitation. As the expansion of the Universe is accelerating and 
general relativity remains unquantized, alternate approaches may prove fruitful. The formalism introduced in this 
paper lends itself quite naturally to this class of theories. 
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APPENDIX A: REVIEW OF GEOMETRIC (CLIFFORD) ALGEBRA 

Many excellent introductions to geometric (or Clifford) algebra exist (for example pi I27I l35l l36|] as well as more 
detailed mathematical treatises (see Q and the various conference proceedings listed in the bibliography), so we will 
make no attempt to be exhaustive. Rather, this appendix and the one following are intended to cover the basics in a 
small amount of space, to get the interested reader started. 

1. Geometric Algebra 

The geometric (or Clifford) product is associative and distributive over addition: 

a(bc) = (ab)c 

a(b + c) = ab + ac (Al) 
and provides a measure of the length of a vector: 

a 2 = aa=\a\ 2 . (A2) 
Thus, if we define a vector a = b + c, using the axioms above, we obtain 

a 2 = (b + c) (b + c) = bb + cc + be + cb = b 2 + c 2 + be + cb. (A3) 
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FIG. 1: Vectors and bivectors illustrating the outer product. 



Since the left side is a scalar, the right side must be as well, so we define the inner product of vectors to be the 
symmetric sum 

b ■ c = ^(bc + cb) = c - b (A4) 

which vanishes when b is orthogonal to c. 

The outer product is then the antisymmetric piece of the geometric product: 

b A c = - (be - cb) = c A b; (A5) 

it is identical to the outer product in differential geometry. Its value is neither a vector nor a scalar, but an oriented 
plane segment called a bivector, whose sides are defined by the vectors making it up (see Fig. ^) . The outer product 
is associative — 

a A (b A c) = (a A 6) A c (A6) 

— and vanishes if the vectors are linearly-dependent. 

Thus, the geometric product of two vectors is the sum of the inner product and outer product — that is, the sum of 
a scalar and a bivector: 

ab = ab + ahb = ba-bha^ ±ba. (A7) 

When we add objects of different types — called grades — together, we have a general object called a multivector. We 
denote a general multivector as the sum of its grades: 

A = (A) + (A) 1 + (A) 2 + ... = Y,(A) 1 (A8) 

7 

where (A) Q is the scalar part, (A} 1 is the vector part, and so on. The grade subscript 7 is an indication of how many 
vectors are needed to construct the multivector via the outer product. 
Therefore, we write 

a b = (ab) and a A b = (ab} 2 (A9) 
for the product of two vectors. We can extend the geometric product to cover arbitrary multivectors, so that 

a (A) =a- <A) r + aA(A) r (A10) 



where 



a.(A> r = i(aA-(-irAo) (All) 



has a grade of r — 1 (and vanishes if r = 0), while 



aA(A) r = i (aA + Aa) (A12) 



has a grade of r + 1. 
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The geometric product also covers multiplication of arbitrary multivectors. We are still most interested in the inner 
and outer products: 

(A) r .(B) s HAB) |r _ s| 

(A) r A (B) s = (AB) r+s , (A13) 

but we also define the commutator product: 

AxB^i(AB-BA) (A14) 

which should not be confused with the vector cross product. It obeys the Jacobi identity 

A x (B x C) + B x (C x A) + C x (A x B) = (A15) 
and is "distributive" over the geometric product: 

A x (BC) = (A x B)C + B(A x C). (A16) 
The commutator product is particularly useful with bivectors, since it preserves the grade of the object it multiplies: 

(B) 2 x(C) r = (BxC) r . (A17) 
One final operation of note is the reverse, which is defined as 

AB = BA (A18) 

for all multivectors, while 

a = a (A19) 

for all vectors. Thus, bivectors change sign under reversion, which geometrically means that their orientation is 
reversed. 

2. Spacetime Algebra 

Up to this point, we have not specified a dimension for our space. Since this paper deals with spacetime, we need 
a space of four dimensions, with a Lorentz signature (i.e., the metric's trace is —2, for consistency with spinorial 
methods). In geometric algebra, this space is denoted by Cl\$ .2], and it is frequently useful to define (locally, at 
least) a set of basis vectors {7^}, which satisfy 

7„ • lu = V»u (A20) 

where r\ — diag(+l, —1, —1, —1) is the metric for a flat Lorentzian spacetime. (Other frames of vectors can be defined; 
a more precise definition of frames is laid out in Section III Dl as well as in Q.) Due to their linear independence, we 
know that 

7«A7 A A7 M A 7l/ A 7c =0 (A21) 
since there do not exist five linearly-independent vectors in the spacetime. In addition, 

7k A 7a A 7 P A 7„ = (A22) 
where e K A^ is the completely antisymmetric Levi-Civita symbol, while 

i = 7o A 7i A 72 A 7 3 = 7o7i7 2 7 3 ( A23 ) 

is the unit tetravector, or pseudoscalar, for spacetime. It has the special properties 

i 2 = -l and i(A) r = (-l) r (A) r i (A24) 

for grade-r multivectors. In addition, the pseudoscalar can be used to define the inner product from the outer product 
1: 

a-B = -aA(Bi)i (A25) 

for any vector a and generic multivector B; the more general version of this relationship demonstrates that the 
existence of a pseudoscalar is equivalent to the existence of a metric. 

The reader may have noticed a strong similarity between the algebra of the basis vectors 7^, and the algebra of 
the Dirac matrices 7^. Indeed, the Dirac matrices are a representation of the Clifford algebra Ci\^, which describes 
spacetime. For more information about matrix representations, please see 0- 
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3. The Even Subalgebra C£^ 3 

The even multivectors (scalars, bivectors, and pseudoscalars) form a subalgebra, which we call the even subalgebra, 
and denote with C£^ 3 . As an algebra, it is isomorphic to the three-dimensional geometric algebra and it is 



possible to perform many spacetime calculations in this restricted space; see |2j, |22j, |24j for more information. 
Proper, orthochronous Lorentz transformations can be written as 0] 

/1(A) = UAU = c B/2 Ac- B/2 (A26) 

where A is an arbitrary multivector, B is a bivector, and the exponential U = exp(B/2) is defined (as in matrix 
algebra) by its Taylor series: 

e B = 1 + B + ^B 2 + ^B 3 + . . . . (A27) 
2 o 

The quantity 

B 2 = B- B + BAB (A28) 

is in general a scalar added to a pseudoscalar, since in four dimensions the exterior product of bivectors cannot be 
assumed to vanish. Thus, this exponential is always a member of the spin group Spin, (1,3) ~ SL(2,C), whose 
Lie algebra is isomorphic to the even subalgebra Ci[ z This subalgebra consists of the scalars, bivectors, and 
pseudoscalars: 

±=(±) Q + (±) 2 + (±) 4 (A29) 

and is closed under the geometric product; all even multivectors commute with the pseudoscalar element as well. We 
call any object in this subalgebra a (Dirac) spinor, following the convention of Hestenes [25l l3fil| . 
The reversion operation allows us to split the algebra into (complex) scalar and bivector pieces: 



(±)b = \(±-±)=(±) 2 -> ( A3 °) 
combining these two operations allows us to decompose any spinor into members of four subspaces: 

t = (t)xs + <£>*B + (±)*B + (ths > (A31) 

where 

(t)m =\{t + t i+ t + i) ( A32 ) 
is the real scalar part, and so forth for all combinations. 

APPENDIX B: GEOMETRIC CALCULUS AND TENSOR ALGEBRA 

The general results of this section are adapted from Q . 

1. Directional and Vector Derivatives 

Consider an arbitrary function of a vector f(a). The rate of change of this function along another vector b is the 
directional derivative, which we write as 

b- d a f{a) = lim -f(a + eb) - f(a); (Bl) 

the notation d a indicates that the derivative is with respect to the vector a. This derivative is scalar-valued, so it 
preserves the grade of the object on which it acts; we also note that when f(a) — a, we obtain b ■ d a a = b. The 
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coordinate derivative is the most useful special case of the directional derivative, and the one we use most often in 
this paper; we find it from letting a = x and b = g^, for some coordinate frame: 

d^(x) = g^d x f(x)=g^df(x). (B2) 

As indicated here, we let 9 = d x . (Equivalently, vectors are frequently defined as directional derivatives.) The 
directional derivative is linear, and obeys the Leibniz rule: 

a-d(A + B) = adA + adB 
(a + b)-dA = a dA + b dA 

a dAB = (a ■ 9A)B + A(a ■ dB). (B3) 

We have implicitly defined d a in this discussion, but its explicit form requires the introduction of some arbitrary 
basis e p , and the use of Eq. (fBlfl : 

8 a = e"e M ■ d a . (B4) 

In turn, this allows us to write the gradient 

d = d a a- d, (B5) 

which we can expand and evaluate using all the previous relations. The operator d a is vector-valued, and therefore 
does not commute with arbitrary multivectors. Due to this, we will find the following identities useful: 

d a ((A) r -a) = (-ir+ 1 9 a (a.(A>„,) = (-ir+V(A) r 
d a ((A) r Aa) = (-l) f fl B (oA(A) r ) = (-in4-r){A) r 

o «A) r a) = d a ((A) r ■ a) + d a ((A) r A a) = (-lf(4 - 2r) <A) r (B6) 

where A is independent of a. Note that these also cover scalars (r = 0), so that 

d a a = d a ■ a = 4 and d a A a = 0. (B7) 



2. Tensors 



Tensors are linear functions of one or more vector or (more rarely) multivector variables. The range of the tensors 
may contain any multivector, but we will most often deal with tensors that preserve the rank of their argument; this 
class covers the vierbein fields, the curvature tensors (Riemann, Weyl, and Ricci), and Lorentz transformations; the 
only real exception we deal with in this paper is the connection bivector, which takes a vector argument and returns a 
bivector. (The connection is usually not considered a tensor, since it transforms inhomogeneously under the Lorentz 
group; however, Sl(b) is linear in its argument, and so is a tensor in that sense. See Appendix|C|) Since tensors are 
linear, they obey the important relation 

T(aA + /3B) = aT(A) + (3T(B), (B8) 

where a and /3 are scalars, and A and B are multivectors. By the same token, for any scalar a, T(a) = a. 
Now let us consider the class of tensors that map vectors onto vectors: 

b' = L(b), 

where we use an underscore to indicate that grade is preserved. The components of the tensor with respect to an 
arbitrary basis {e^} are 

V = e„ ■ T{e v ), (B9) 

so that 

T{a) = a»T(e^ = o"T\e A . 
The adjoint (or transpose) of a tensor is found from the definition 

a ■ 7(6) = b ■ T(a) -> 7(6) = d a a ■ 7(6), (BIO) 
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where d a is defined as in the previous section. A tensor of this type can be extended to operate on a general multivector 
via the outermorphism: 

n n 

where 

n 

A aj ■ = a : A a 2 A . . . A a n . (B12) 

i=i 

(Note that curvature, although it takes a bivector argument, is not an outermorphism: it cannot be decomposed as 
R(a A 6) = R'(a) A R'(b).) Since the pseudoscalar i is unique to a multiplicative scalar, the outermorphism acting on 
a pseudoscalar must be a defined quantity: 

T(i)=idetT -> det T = -iT(i), (B13) 

where det T is the determinant of the tensor. If the determinant is nonzero, we can write the inverse of the operator 
as 

T _1 (a) =T(oi)(idetT)- 1 . (B14) 

The trace of a tensor of this type is defined as 

trT = a ■ 7(a) = e" • 7(e„) = e" • e A 7% = T%, (B15) 

where we have defined a basis for evaluation. When operating on a tensor of a more general type, this is called the 
contraction: 

d a -T(a,b,...) = S(b,...). (B16) 
The protraction and traction are defined analogously, using the outer and geometric products, respectively: 

9 a A T(a, &,...) d a T(a,b,...) =d a - T(a,b, . . .) + d a A T(a,b,...). (B17) 
The main place where these operations are used in this paper is in defining the Ricci tensor, 

R(a) = d b R{b A a) = d b ■ R(b A a) -> d b A R{b A a) = 0, (B18) 
and in noting that the Weyl tensor is "tractionless" : 

d b C(bAa)=0. (B19) 

APPENDIX C: DERIVATION OF THE CONNECTION 

As discussed in Section III Dl a coordinate frame satisfies Eq. (|2.4U|) : 

VAfif" = -> V Ag» = -d a A(n{a)-g»). (CI) 

In geometric algebra, it is possible to solve Eq. H(J1|) for the connection bivectors; we are not aware of a solution in 
any other language. 

To isolate fl(a), we first "differentiate" out g^: 

9„ A V A g» = -g M A d a A (n(o) • g^) = d a A ( 9fi A (fl(o) • g^)) = -29a A O(o). (C2) 
After noting that 

= V A (g v A g v ) = g^ A V A - A V A 9fl , (C3) 
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we proceed by taking the inner product of the left and right sides of Eq. ljC2fl with an arbitrary vector 6 
b-tAVAg") = 6- (>AVAgJ = 6"V A - A (b-Vg.) +g"A S A 6- fv A g 



= ^VAg„-g"A [b-Vg^ +g"AV&, + VAfc (C4a) 
b ■ [d a A 12(a)] = 12(6) + c? a A (12(a) ■ b) = 12(b) + VA6-VA6. (C4b) 

(For simplicity, we have introduced the shorthand Va = g A • V, 6 M = b ■ g M and = b ■ g .) We can evaluate the 
covariant curl of b by considering its components b = b u g v : 

V A b = g x A V A (6„ff ") = g x A { (VaMa " + b v Vxg" } = -g v A V6„ (C5) 

due to the torsion-free condition and the fact that the components are scalars. 
Solving Eq. fH4)) for 12(6) yields 

212(6) - V A 6 + g" A V6 M - fe^V A g M + g" A (ft • Vg J . (C6) 
This expression is further simplified by unwrapping the curl: 

V A 6 = g x A Va (6" ff/t ) = b»V Ag^-g^A V6". (C7) 
When plugged into Eq. I|C6(I . all of this gives us our final expression: 

«(&) = ^(/A%-^AW+/A(b.V 9fI 

= ^ {g M A g A 6 ■ (v A gJ - g M A g A 6 • (v Aff ") + g» A (b ■ Vg J } . (C8) 

The second form of this demonstrates the linearity in 6, while the first is useful in the case 6 = g a , where the 
components are 6 M = 5£ and 6 M = g^ a ; here, Eq. I|C8|I reduces to Eq. I|2.41bj) : 

12 Q = 12(g Q ) = \{g v A Vg m + g M A V Q g J . (C9) 

This last equation is derived in Snygg using Christoffel symbols and Ricci coefficients, while the more general 
form is used in Doran |31| without referring directly to a coordinate frame. 

APPENDIX D: NULL ROTATIONS 

In a physical problem, we may have an obvious null direction specified; to work out the symmetries of the problem, 
it would be useful to perform an operation that preserves that direction, but changes all the others. A null Lorentz 
transformation is the operation we are looking for, and is defined in C£i t 3 by 

/1(A) = RAR (Dl) 

where 

R = e f/2 = 1 + -f , (D2) 

and f 2 = is a bivector. To preserve the null vector I — 2lt (and hence I), 

f = p A I where p ■ I = (D3) 

and p is an undetermined spacelike vector. Due to the orthogonality condition, only two components of this vector 
need to be determined, which we will write as the complex component of the bivector f : 

f = V2aZ = 2(pt) B (D4) 
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so that a = p u — ip v in terms of the bivector basis {e, u. v}. Under the null rotation, we have a new set of tetrads 

V = l 

c^rn) + aocl 

_ aU. (D5) 

Under this same transformation, the bivector basis becomes 



Z' = -V2(iW> =Z 



ii 



z'i = (L'm') B - (m'm'^B = Zi + V2aZ 

Z 2 = -V2(n'rn,' 1 ) B = Z 2 - V2aZi - a 2 Z a . (D6) 

In this paper, we will not need the complementary transformation, in which n is preserved; however, it is simple to 
obtain those relations using an analogous discussion |18| . 

Under this transformation, the Weyl matrix elements as defined in Eq. (|6.16l) change as follows: 

(D7a) 
(D7b) 
(D7c) 
(D7d) 
(D7e) 

where the Z' a are given in Eq. l|D6|l and a is the complex scalar defined in Eq. I)D4|) . Using transformations like this 
(and ones that preserve n), we can adapt our system of tetrads to seek out the Petrov classification for an arbitrary 
geometry. As noted in [lg, we can also find the classification of a known tensor by solving the equation ty' 4 = 0, 
which is an algebraic equation in a. 



% = 


-{Z' C(Z' )) S 


= *0 




*i = 


-(Z' C(Z[)) S 


= *i + V2a^ 






-(zic(z;)) s 


+ (z' Q c(zy 2 ) s = 


^2 + iV2a^i + 3a 2 * 


*3 = 


-(z;c(z 2 )) 5 


= *3 - V2a^> 2 - 


- 3a 2 *i - V2a 3 * 


% = 


-<Z_C(Z' 2 )) S 


= * 4 - 2\/2a^ 3 
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